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An  expression  for  the  emitted  fluorescence  radiation  by  an  atom  near  a 
phase  conjugator  is  derived  from  a  plane-wave  expansion  of  the  electric  field 
operator  and  Heisenberg's  equation  of  motion  for  the  annihilation  operator. 
The  result  is  compared  to  a  solution  which  was  found  previously,  based  on  the 
classical  Maxwell  equations.  It  is  shown  that  both  theories  yield  the  same 
expression  for  the  field  in  the  radiation  zone,  in  the  limit  of  a  transparent 
medium.  This  confirms  the  correctness  of  either  approach  to  the  problem  of 
optical  phase  conjugation  of  atomic  radiation. 
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I.  INTRODUCTION 

In  a  previous  paper*  we  have  studied  the  behavior  of  an  atom  near  the 
surface  of  a  four-wave  mixing  phase  conjugator  (PC).  We  calculated  the 
fluorescent  photon  emission  rate  by  solving  Maxwell's  equations  for  a  dipole 
H  near  a  PC.  The  electric  field  E(r,t)  was  interpreted  as  a  quantum  operator 
field,  with  the  argument  that  the  classical  Maxwell  equations  must  be 
identical  in  form  to  the  (quantum)  Heisenberg  equations  of  motion  for  the 
operator  fields.  Nevertheless,  questions  can  be  raised  as  to  whether  the 
identification  of  Maxwell's  equations  with  Heisenberg's  equations  can  be 
justified  for  the  problem  under  consideration.  In  particular,  the  pump  beams 
(with  frequency  u)  for  the  four-wave  mixing  process  are  taken  into  account 
parametrically,  and  are  represented  by  classical  plane  waves.  This  leads  to 
factors  of  exp(-2iwt)  in  the  expression  for  the  fluorescence  radiation  field, 
rather  than  a(t)  ,  with  a(t)  the  annihilation  operator  for  a  photon  in  the 

pump  beam.  In  addition,  our  results  are  in  conflict  with  the  results  of 

2  3  ( - ) 

Hendriks  and  Nienhuis,  ’  who  did  not  find  the  terms  proportional  to  c ) 

(raising  part)  in  the  expression  for  the  fluorescence  field. 

In  this  paper  we  solve  the  Heisenberg  equation  for  the  annihilation 
operator  a^(t)  for  a  photon  with  wave  vector  k  and  polarization  a  (either  s  - 
surface  polarized  or  p  -  plane  polarized).  The  solution  is  applied  to 
evaluate  the  fluorescent  radiation  field  int  he  far  zone.  These  calculations 
are  completely  independent  from  our  previous  method,  and  the  results  can  be 
used  to  verify  the  consistency  of  our  approach. 

II.  ELECTRIC  FIELD 

The  e lec tr ic - f ie Id  operator  E(r,t)  can  be  represented  as  a  sum  of 
polarized  plane  waves.  In  the  region  z  >  0  (above  the  PC,  where  the  atom  is), 


we  have  incident  (inc)  waves  with  wave  vector  k  and  polarization  a.  These 
waves  gives  rise  to  specular  ly  -  ref  lected  (r)  waves,  and  to  phase-conjugated 
(pc)  waves  which  travel  in  the  direction  of  -k.  In  addition,  there  are 
transmitted  (t)  waves  which  have  their  origin  in  waves  which  are  incident  on 
the  medium  from  the  other  side  of  the  PC.  The  four-wave  mixer  produces  also  a 
nonlinear  (ni)  wave  in  z  >  0. 

The  amplitude  of  each  of  the  generated  waves  is  related  to  the  amplitude 
of  the  corresponding  incident  wave  by  a  Fresnel  coefficient.  The  electric 
field  is  explicitly, 
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(2.1) 


where  a  prime  (double  prime)  on  the  summation  sign  indicates  a  sura  over  waves 

which  propagate  in  the  z  <  0  (z  >  0)  direction  only.  The  unit  polarization 

vectors  e^  ^  and  the  wave  vectors  k.  are  defined  in  Appendix  A.  The 
-  i  ~ l 

Heisenberg  operator  ako(t)  ls  the  annihilation  operator  for  a  photon  which 


is  incident  on  the  PC,  either  from  z  >  0  or  from  the  other  side  (back  port)  of 


the  medium.  We  shall  take  the  Schrodinger  picture  and  Heisenberg 
coincide  at  t  -  0. 


III.  EQUATION  OF  MOTION 

The  only  unknown  in  expression  (2.1)  for  E(r,t)  is  the 
operator  a^(t)  for  t  >  0.  Its  equation  of  motion  is 


'[ a,_  (t),  H] 


with  initial  condition  a.  (0)  —  a.  The  Hamiltonian  H  can  be  wri 
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for  the  Hamiltonian  of  the  radiation  field,  and  where  w  -  ck. 

Hamiltonian  H  can  remain  unspecified  because  it  commutes  with 
a  r 

the  interaction,  we  take 
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The  atomic 

a.  (t).  For 
ka 


(3.4) 


in  terms  of  the  atomic  dipole  moment  operator  m  and  the  position  h  -  he^  of 
the  atom.  With  Eqs .  (3.3)  and  (3.4)  the  equation  of  motion  becomes 


itf  ^  a^ff(t)  -  a^ff(t)  -  M(t)» [a^ff(t) ,E(h, t) 


(3.5) 


where  E(h,t)  follows  from  Eq .  (2.1)  with  r  —  h. 


V.  SOLUTION 


An  integral  of  Eq.  (3.5)  is 
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.  (4.1) 


When  we  subsitute  Eq .  (2.1)  for  E(h.t'),  then  it  appears  that  many  terms  in 
the  integrand  oscillate  at  optical  frequencies,  and  these  terms  average  out  to 
zero  on  a  time  scale  of  an  optical  cycle.  The  time  dependence  of  the 
annihilation  operator  would  be  ak^(t)  -  a^exp(-iwt)  in  free  evolution,  which 
contains  only  a  single  positive  frequency.  When  the  interaction  is  taken  into 
account,  a^(t)  acquires  a  spectral  width  around  the  frequency  w,  but  it  is 
still  a  pos i t ive - f requency  operator.  Similarly,  when  we  split  the  dipole 
moment  into  a  positive  and  a  negative  frequency  part  as 


M(t)  -  M(t)(+)  + 


(4.2) 


then  ^(t)^'*  -  r*(  +  ^exp(-iw  t)  in  free  evolution  (for  a  two-level  atom  with 

transition  frequency  w  ).  The  third  time  dependence  enters  as  exp(±2iwt). 

o 

Then  we  drop  all  terms  in  the  integrand  which  oscillate  with  optical 


o 


rrequencies,  and  retain  the  terms  which  oscillate  with  the  difference  of  two 
optical  frequencies.  Furthermore,  we  notice  chat  the  first  term  on  the  right- 
hand  side  of  Eq .  (4.1)  yields  the  vacuum  field  E  (r,t). 

Therefore,  we  can  write 


E(£.t)  -  Ev(r,t)  +  E  (r ,  t) 


(4.3) 


where  E^Cr.t)  is  given  by  Eq.  (2.1)  with  a^  (t)  -*■  a^exp  ( -  iwt )  .  The  source 
rield  E  (r, t)  follows  from  the  second  term  on  the  right-hand  side  of  Eq. 
4.1).  Ve  find  explicitly, 
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in  terms  of  the  operators 
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V.  ASYMPTOTIC  EXPANSION 

In  order  to  simplify  the  solution  (4.4)  for  Eg(r,t)  we  consider  its 

value  in  the  radiation  zone  (r  -  <*>  and  z  >  0).  The  operators  bk^(t)  and 

bk^(t)  are  independent  of  r,  so  that  all  r-dependent  factors  appear  in  the 

form  exp(ik^*r).  These  exponentials  are  multiplied  by  a  function  of  k^,  and 

the  result  is  summed  over  all  values  of  k  .  Therefore,  all  terms  in  Eq.  (4.4) 

-a 

have  the  generic  form  V  3  ^  g(k^)exp( ik^> r) .  The  summation  runs  either  over 
wave  vectors  with  only  positive  z-components  or  over  wave  vectors  with  only 
negative  z-components.  Changing  the  summation  into  an  integration  gives 
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where  the  superscript  (±)  on  the  region  of  solid  angle  indicates 

the  z-components  of  the  wave  vectors  k  .  Then  we  can  make 

~a 

expansion  of  the  angular  integral  with  the  methed  of  stat io».cu.y 
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result  is 
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'-'here  e  is  the  radial  spherical  unit  vector  which  points  in  the  observation 
-r 

direction.  This  direction  will  be  specified  by  the  spherical  angles  6  and  <j>  . 
The  asymptotic  expansion  effectively  filters  out  the  value  of  g(k^)  for  k^  -  ± 
k  e  corresponding  to  the  plane  wave  exp(ik?*r)  which  travels  into  the 
observation  direction. 

with  Eqs.  (5.1)  and  (5.2),  the  asymptotic  expansion  of  E  r,t)  is  round 
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’I .  POLARIZATION  VECTORS 
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exDressed  in  the  sDherical  unit  vectors  e,  and  e,.  These  calculations  are 

-4  -4> 

similar  'but  not  identical)  to  the  corresponding  calculations  in  Ref  1  (Sec. 
IX) .  Here  it  is  convenient  to  introduce  the  notation 


e  -  +e  . 

-S 


(6.1) 


■hich  will  enable  us  to  express  the  results  for  s-waves  and  p-waves  m  a 
;  ir.sie  formula.  Xe  find 
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and  similar  expressions  for  X  .  X  and  X  .  Apparently,  the  term 

pc, a  t ,  a  nr ,o 

orooortional  to  X  .  is  an  advanced  (non-causal)  contribution  to  E  , 

adv, a  ~s 

corresponding  to  incoming  spherical  waves.  The  other  four  terms  are  retarded 


:ausal)  solutions,  and  they  have  the  form  of  outgoing  spherical  waves. 
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Next  we  substitute  expressions  (4.5)  and  (4.6)  for  b^  (t)  anc*  b/^Ct), 

re soec tive Ly ,  into  the  results  for  the  X  's.  Then  we  evaluate  the 

ct  ,  a 

polarization  vectors  in  b,  (t)  and  b,1  (t)  at  the  indicated  wave  vectors. 

r£  O  K.  G 

After  lengthy  calculations  we  find 
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and  four  other  similar  expressions.  Here  we  introduced  the  abbreviations 
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and  ail  Fresnel  coeff icients ,  like  P  ,  R  ,  etc.,  are  evaluated  at  frequency  u 

a  a 

and  an  angle  of  incidence  equal  to  the  polar  angle  ( 8 )  of  observation. 


VII.  TVO- LEVEL  ATOM 

So  far  we  have  not  made  any  assumption  about  the  atom  and  its  dipole 
moment  n(t) .  In  order  to  evaluate  the  expression  for  E  (r,t)  further,  we 
assume  that  only  two  atomic  levels,  which  might  be  degenerate,  are  of 


11 


relevance.  The  level  separation  is  •  In  free  evolution,  the  positive  and 
negative  frequency  part  of  n(t)  obey  the  identity 
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We  shall  use  this  as  an  approximation  in  equations  of  the  type  (6.5).  Then, 
for  instance,  Eq.  (6.5)  can  be  written  as 
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In  every  term,  both  the  exponential  in  the  integrand  and  the  exponential  in 
front  of  the  integral  have  the  same  frequency.  When  the  distance  r  is  much 
larger  than  a  wavelength  and  the  retardation  r/c  is  much  larger  than  an 
optical  cycle,  then  the  integrals  can  be  approximated  by^ 


- iw  (t±r-r/c)  pt  iu  t' 
e  a  dt'  e  a  -2k  6 (w  ) 

Jo  a 


(7.3) 


In  writing  out  all  ^  as  in  Eq .  (7.2),  it  can  be  shown  that  every 

integral  can  be  written  as  in  Ea.  (7.3).  The  only  exception  is  X  ,  ,  which 

adv ,  a 

contains  the  integral  of  Eq.  (7.3)  with  r  replaced  by  -r.  This  makes  the 

right-hand  side  zero,  instead  of  2n6(u  ),  as  shown  in  Ref.  5.  Therefore, 

a 


X  ,  -  0 

adv,  a 


(7.4) 


as  it  should.  Then  we  set  k  -  w/c  and  carry  out  the  integrations  over  w. 
This  gives,  for  instance, 


X 

?  ,a 


2  rru" 


[R  e+ • n( t- t -r/c) ^ 

C7— <7  — 


|R  I  ^e+>n(t+T -r/c) ^  +  ^ 
a 1  ~a  ~ 


— 

+  R  P 


aPff  ^^(t-r-r/c)(-)e-2iw(t-r-r/c) 


(7.5) 


where 


P 

a 


P  (w  ,s) 
a  o 


P  -  P  (2u-<w  ,  8 ) 
o  a  o 


(7.6) 

(7.7) 


and  similarly  for  other  Fresnel  coefficients. 

VIII.  TOTAL  SOURCE  FTF.T.D 

The  total  field  Eg(r,t)  has  four  contributions  of  the  form  (7.5), 
according  to  Eq.  (6.2).  We  set  t  -  r/c  -*  t,  and  introduce  the  polarization 
vectors  (without  the  ±  superscripts) 
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e  -  e  , 
-s  ~<j> 


e  -  e 
-p  ~d 


(8.1) 


€  -  e.  «  -e,  -  2sin0e 

-s  -4>  ~p  -3  -2 


(8.2) 


rhen  we  define  the  parameters 


|R  I2  +  t T' 1 2 
‘  cr 1  1  a  * 


I?/  •  !»;l2 


(8.3) 


_  *  _  _  *  * 
R  P  -  R  P  +  T'N'  -  T'N' 
a  a  a  a  a  a  a  a 


(8.4) 


which  will  depend  on  the  angle  d,  in  general.  When  we  group  together  all 
terms  with  and  all  terms  with  p(t)^  ^ ,  then  the  field  assumes  the 
remarkably  simple  form 


2 

cj  e 


Es(r,t) 


5 - - —  \  e  (a+.p(t)(  +  )  +  e-2iwta‘^(t)('))  +  H.c.  .(8.5) 

2  /  —a  -a  -  ~ 

itf  c  r  L- 


o  a 


in  terms  of  the  polarization- like  vectors 


2iw  T 

+  a  o  „ 

a  -  y  e  +  e  Re 

-a  a~a  a -a 


(8.6) 


,  2i«  r  . 

*  0  1 
a  —  -Pe+e  ye 
-o  a-o  0-0 


(8.7) 


Alternatively,  we  can  group  the  terms  as 
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w2e'1Wor 

Es(r,t)  -  -21-—  KM, •£,)£,  +  (M  -e  )e  >  +  H.c. 

4fff  c  r 

o 


(3.8) 


with 


M  -  7%<t)(  +  )  +  e  °  R  /i'(t)(  +  ) 
-®  P~  P~ 


. . —  ,  .  .  2iwr.  .  . 

-2lwt.  *  ,_.(-)  o  l  .(-K 

-  e  (P  m(C)  -  e  7  M  (t)  ) 

p-  p- 


(8.9) 


.  ,  .  2iw  r  ,  . 

M  -  7aM(t)(+)  -  e  °  R  M' (t)(+) 
-0  s—  s  — 


2iu  r  . 


e-2iwc(p*^(t)(-)  +  e  o  (t)^*)) 


(8.10) 


The  mirror  dipole  n'  is  defined  as 


A*'  (t)  -  Mj.(t)  '  if ||  (c) 


(8.11) 


in  terms  of  the  perpendicular  and  parallel  components  of  ff(t)  with  respect  to 
the  surface  z  -  0.  Expressions  (8.8)  -  (8.10)  would  be  identical  to  our 
earlier  results  (12.5)  -  (12.7)  from  Ref.l.  if  the  parameters  7^  and  7  would 


be  ya  -  1  and  yL  -  0. 
a  a 


We  have  not  used  any  of  the  properties  of  the  Fresnel  coefficients  in 
the  derivation  of  the  results  of  the  previous  section.  These  coefficients  can 
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be  calculated  explicitly, '  but  the  result  is  very  complicated.  Therefore,  we 
consider  two  limiting  cases  of  practical  interest. 

A.  Dielectric  Laver 

When  we  turn  off  the  pump  beams,  then  the  medium  becomes  an  ordinary 
dielectric.  Therefore, 


P  -  N  -  0 
a  a 


which  gives 


The  non-zero  Fresnel  coefficients  are  related  by 


so  that 


(9.1) 


(9.2) 


(9.3) 


This  shows  that  for  a  dielectric  layer,  the  results  (8.8)  -  (8.10)  are 
identical  to  those  in  Ref.l.  Notice  that  the  term  which  is  proportional  to 
^(t)<‘  ^  in  Eq.  (8.5)  disappears  in  this  limit. 
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B .  Transparent  PC 

When  the  dielectric  constant  equals  unity,  the  specular  waves  vanish. 
This  gives 


R  -  N  -  0  (9.5) 

a  a 


The  non-zero  Fresnel  coefficients  are  now  related  by 


(9.6) 


which  holds  for  any  polarization,  angle  of  incidence  and  frequency, 
Therefore,  we  find  again 


7*  -  1  ,  7^-0  . 


(9.7) 


X.  CONCLUSION 

We  have  derived  an  expression  for  the  fluorescence  radiation  field  which 

is  emitted  by  an  atomic  dipole  near  the  surface  of  a  PC.  The  starting  point 

was  the  standard  plane-wave  expansion  of  the  electric  field  in  terms  of 

annihilation  and  creation  operators.  Then  we  solved  the  Heisenberg  equation 

of  motion  for  the  annihilation  operator.  This  gives  rise  to  two  contributions 

to  the  electric  field:  the  vacuum  field  E  and  the  source  field  E  ,  which  is 

~v  -s 

generated  by  the  dipole.  The  form  of  E^  follows  trivially  from  the  choice  of 
H^,  but  the  form  of  E^  depends  on  the  choice  of  interaction  Hamiltonian  and 
the  structure  of  E^  (which  equals  E(r,0)).  We  have  shown  that  Eg  is  identical 
to  the  solution  of  Maxwell's  equations,  as  found  previously,  provided  that 
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a 


1 


(10.1) 


These  parameters  do  not  appear  in  the  solution  of  Maxwell's  equations,  which 
indicates  that  both  approaches  are  independent  indeed.  The  parameters  7a  and 

7“  are  determined  bv  the  Fresnel  reflection  and  transmission  coefficients  for 

a  - 

a  plane  wave,  and  they  depend  on  the  polarization,  frequency  and  angle  of 
incidence.  The  general  form  of  the  Fresnel  coefficients  is  extremely 
complicated,  which  prohibits  the  verification  of  Eq .  (10.1)  for  the  most 

general  case.  We  have  shown,  however,  that  for  a  dielectric  layer  and  for  a 
transparent  medium  the  relations  in  Eq .  (10.1)  are  satisfied,  which  covers 

most  cases  of  practical  interest. 
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APPENDIX 

Here  we  summarize  che  expressions  for  the  various  wave  vectors  in  Eq. 
(2.1)  and  the  phase  conventions  for  the  unit  polarization  vectors.  Any  wave 
can  be  decomposed  as 


k.  -  ku  +  k.  e  ,  (A. 1) 

-i  -||  i,2  -z 

and  any  k^  in  Eq .  (2.1)  which  corresponds  to  a  given  incident  wave  must  have 
che  same  parallel  component  k.j  with  respect  to  the  surface  of  the  medium.  The 
polarization  vectors  for  s-waves  and  p-waves  are  chosen  as 


x  e 
-z 


(A. 2) 


x  e, 
-k^s 


(A. 3) 


respectively.  An  incident  wave  from  z  >  0,  and  with  wave  vector  k  -  k||  + 
k^e  generates  a  specular  (r)  wave  and  a  phase -conjugated  (pc)  wave.  The  z- 
components  of  their  wave  vectors  are  determined  by  the  dispersion  relation, 
and  found  to  be 


(A. 4) 

(A. 5) 

in  terms  of 


k  -  -  k 
r  ,z  z 


,  PTl  72 

k  -  -  vp  k  -  kii 

pc.z  F  II 


p  -  (2w  -  w)/w 


(A. 6) 
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A  wave  with  wave  vector  k  which  is  incident  on  the  layer  from  z  <  0,  generates 
a  transmitted  (t)  wave  in  z  >  0  with  wave  vector  k.  In  addition,  a  nonlinear 
(ni)  wave  is  produced  in  z  >  0,  which  has  the  same  frequency  shift  with 
respect  to  the  incident  wave  as  the  pc-wave. 
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